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Global Formulations of Lagrangian and Hamiltonian Mechanics 

on Two-Spheres 

Taeyoung Lee*, Melvin Lcok', and N. Harris McClamroch 


Abstract —This paper provides global formulations of La¬ 
grangian and Hamiltonian variational dynamics evolving on 
the product of an arbitrary number of two-spheres. Four types 
of Euler-Lagrange equations and Hamilton’s equations are 
developed in a coordinate-free fashion on two-spheres, without 
relying on local parameterizations that may lead to singularities 
and cumbersome equations of motion. The proposed intrinsic 
formulations of Lagrangian and Hamiltonian dynamics are 
novel in that they incorporate the geometry of two-spheres, 
resulting in equations of motion that are expressed compactly, 
and they are useful in analysis and computation of the global 
dynamics. 

I. Introduction 

The two-sphere is the two-dimensional manifold that is 
composed of unit-vectors in M 3 . There are a wide variety of 
dynamical systems that evolve on multiple copies of two- 
spheres. In robotics, the configuration of articulated robotic 
arms interconnected by spherical joints is represented by 
two-spheres [1], Since the surface of the Earth is approx¬ 
imately a sphere, spherical dynamics arise readily in earth 
science and meteorology [2]. In quantum mechanics, the pure 
state space of a two-level quantum mechanical system is a 
two-sphere, referred to as the Bloch sphere [3], 

In most of the existing literature on dynamical systems 
evolving on two-spheres, the unit-sphere is parameterized by 
two angles. For example, a point on the two-sphere is often 
described by its longitude and latitude. Parametrizing the 
two-sphere by two angles is straightforward, and the angles 
are typically viewed as being in an open subset of R 2 . 

However, such parameterizations of the two-sphere suffer 
from the following two main issues. First, parameterizations 
represent the two-sphere only locally. This can be easily 
observed from the fact that the longitude is not well defined 
at the north pole and the south pole. This causes a singularity 
in representing the kinematics on the two-sphere, especially 
when converting the velocity of a curve on the two-sphere 
into the time-derivatives of the longitude and the latitude. 
This yields numerical ill-conditioning in the vicinity of those 
singularities, which cannot be avoided unless one switches 
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coordinate charts, which becomes problematic when trying 
to track motions with large angular deviations. 

The second issue is that the equations of motion of 
dynamical systems on the two-sphere become exceedingly 
complicated when expressed using local coordinates and 
necessarily involve complicated trigonometric expressions. 
For example, the dynamics of a multiple spherical pendulum, 
written in terms of angles, is extremely complicated. 

This paper aims to provide global formulations of dynam¬ 
ics evolving on the product of an arbitrary number of two- 
spheres. In particular, we study dynamical systems that can 
be viewed as Lagrangian systems or as Hamiltonian systems 
that encompass a large class of mechanical systems that 
appear in robotics, structural dynamics, quantum mechanics 
or meteorology. Most importantly, the unit-vectors on the 
two-sphere are regarded as elements of a manifold, and 
dynamics are formulated directly on two-spheres in a global 
fashion via variational principles. 

This geometric formulation is said to be coordinate-free, 
as it does not require the use of local charts, coordinates or 
parameters that may lead to singularities or ambiguities in 
the representation. As such, it can be applied to arbitrarily 
large angle rotational maneuvers on the two-spheres globally. 
Furthermore, this provides an efficient and elegant way to 
formulate, analyze, and compute the dynamics and their 
temporal evolutions. The corresponding mathematical model 
developed on two-spheres is nicely structured and elegant. 
This representational efficiency has a substantial practical ad¬ 
vantage compared with local coordinates for many complex 
dynamical systems; this fact has not been appreciated by the 
applied scientific and engineering communities. 

In short, the main contribution of this paper is provid¬ 
ing geometric formulations of the equations of motion for 
Lagrangian and Hamiltonian systems that evolve on two- 
spheres using variational methods. The proposed global 
formulations, that do not require local charts, have not 
been previously studied, even in the well-known literature 
on geometric mechanics, such as [4], [5], [6]. Preliminary 
results have been given in [7], [8], where Euler-Lagrange 
equations are developed for a certain class of mechanical sys¬ 
tems whose kinetic energy is repressed as a quadratic form 
with fixed inertia elements. This paper provides both Euler- 
Lagrange equations and Hamilton’s equations for arbitrary 
mechanical systems without such restrictions. 


II. The Two-Sphere 

The two-sphere is the two-dimensional manifold of unit- 
vectors in R 3 , i.e., 

S 2 = {gGR 3 |||g|| = l}. (1) 


It is composed of the set of points that have the unit distance 
from the origin in R 3 . The tangent space of the two-sphere 
at q G S 2 corresponds to the two-dimensional plane that is 
tangent to the sphere at the point g, and it is given by 

TqS 2 = u G R 3 1 q ■ Z = 0}. (2) 

Throughout this paper, the standard dot product between two 
vectors is denoted by x ■ y = x T y for any x, y G R n . 

Therefore, for any curve q(t) : M —> S 2 on the two- 
sphere parameterized by time t, its time derivative satisfies 
q(t) ■ q{t) = 0. From now on, we do not explicitly denote 
dependence on time for brevity, unless needed. This implies 
that there exists an angular velocity u> : R —>• R 3 such that 


q = oj x q = S(cv)q, 


(3) 


where the hat map S(-) : R 3 —> R 3x3 is defined such that 
S(x)y = x x y and S T (x ) = —S(x) for any x,y G R 3 . 
More explicitly. 


S(u>) 


0 

-w 3 

w 2 

w 3 

0 

-U>1 

—w 2 

Wl 

0 


(4) 


Without loss of generality, the angular velocity is constrained 
to be orthogonal to q, i.e., uj ■ q = 0. Therefore, the three 
vectors q, q, and ui are mutually orthogonal, and the angular 
velocity can be written as 


w = S{q)q. (5) 

It follows that 6j = S(q)q and uj is perpendicular to q as 
well. 


III. Lagrangian Mechanics on Two-Spheres 


Let qi : -A- S 2 be a differentiable curve for i G 

{1,..., n}. The variation is a parameterized curve defined as 

q\ : (—c, c) x [to,tf] —> S 2 for c > 0, such that q°(t) = qt(t) 
for any t G [t 0 ,tf] and gf(f 0 ) = qi(t 0 ), q%(t f ) = q(t f ) for 
any e G (—c, c). 

We can express the variation of the curve qi using the 
matrix exponential map as follows: 

Qi(t) = e £S( ^ (t)) qi (t), ( 6 ) 


for differentiable curves 7 * : [to,tf] —>• R 3 satisfying 
7 i(fo) = 7 i(tf) = 0. Since the exponent eS( 7 *) is skew- 
symmetric, the exponential matrix is in the special orthog¬ 
onal group, SO(3) = {i? G R 3x3 | R T R = J, det[i?] = 1}, 
thereby guaranteeing that the variation is a parameterized 
curve on S 2 , i.e., ||g|(f)|| = 1. There is no loss of generality 
in requiring that 7 j(t) ■ qi(t) = 0 for all to < t < tf, that is, 
7 i and qi are orthogonal. In short, the variation of the curve 
qi in S 2 is expressed in terms of a curve in R 3 via ©. 

The corresponding infinitesimal variations are given by 

Sqi(t) = — q e {t) = 5(7i(f))gi(f), (7) 

de 6=0 

and satisfy Sqi(to) = Sqi(tf) = 0. Since the variation and the 
differentiation commute, the expression for the infinitesimal 
variations of the time derivatives are given by 


5qi(t) 


d 

de 


q e (t) = S('j i (t))q z (t) + S'( 7 i(f))g l (f). 

6=0 


( 8 ) 


These expressions define the infinitesimal variations for a 
vector function q = (gi,...,g n ) : [toRf] —> (S 2 )". The 
infinitesimal variations are important ingredients to derive 
the Euler-Lagrange equations on (S 2 )™. We subsequently 
suppress the time argument, thereby simplifying the notation. 

The action integral is the integral of the Lagrangian 
function along a motion of the system over a fixed time 
period, i.e., © = jf* L(q,q)dt. The infinitesimal variation 
of the action integral is given by 


We consider dynamical systems evolving on the product of 
n copies of two-spheres, namely (S 2 )”. The corresponding 
Euler-Lagrange equations are first obtained in terms of a 
Lagrangian expressed in terms of the configurations and 
the time derivatives of the configurations, namely (g, g). 
A second form of Euler-Lagrange equations is obtained in 
terms of a modified Lagrangian expressed in terms of the 
configurations and the angular velocities, (q,ui). In each 
case, these Euler-Lagrange equations are simplified for the 
important case that the kinetic energy function is a quadratic 
function of the time derivatives of the configurations. 

A. Euler-Lagrange equations in terms of (g, g) 

Suppose that a Lagrangian L(q,q) : T(S 2 ) rl —> R 1 
is given on the tangent bundle of (S 2 )", where (g, g) = 
((gi,..., q n ), (gi,..., q n )) G T(S 2 ) n . For example, it can be 
defined as the difference between a kinetic energy and a po¬ 
tential energy. We derive the corresponding Euler-Lagrange 
equations according to Hamilton’s variational principle. 


d© = 


de 


l s 


L(q e , g £ ) dt = 0. 

e—0 J to 

This can be rewritten in terms of dqi as 

5 © = [ tf -6q t +ir dL{q ’ q) 


'to 


. 2=1 


2=1 


dqi 



dt. 

(9) 


Let fi G M 3 be the generalized force acting on g^ The 
corresponding virtual work is given by 

n 

5W = Y J h-^- ( 10 ) 

i—l 

According to the Lagrange-d’Alembert principle, we have 
<5© = — 6W for any variations. We now substitute the 
expressions for the infinitesimal variations of the motion ([7]i 
and © into this, and we simplify the result to obtain the 
Euler-Lagrange equations expressed in terms of (g, g). 









Proposition 1: Consider a Lagrangian L(q,q) : 
T(S 2 ) n —>• R for a mechanical system evolving on 
the product of two-spheres, with the generalized force /,. 
The Euler-Lagrange equations are given by 

T \ f d (dL{q,q)\ 8L(q,q ) An 

qiq '>\dtX dq t ) % fi j~ ' 

(ID 

for i £ {1.... . n}. Here, the 3x3 identity matrix is denoted 
by I 3 £ R 3x3 . 

Proof: See Appendix IBl ■ 

This describes the evolution of the dynamical flow ( q , q) £ 
T (S 2 )'"' on the tangent bundle of the configuration man¬ 
ifold (S 2 ) n . In the above equation (fill , the expression 
in the braces corresponds to the Euler-Lagrange equation 
for dynamical systems evolving on R". It is interesting to 
note that the Euler-Lagrange equations on S 2 corresponds 
to its orthogonal projection onto the plane normal to qi , 
represented by the matrix /•-> - q,qj. 

Next, we consider the important case that the kinetic 
energy is a quadratic function of q, and the potential energy 
U is given as a function of q, i.e., the Lagrangian can be 
written as 

1 n 

L(q,q) = 9 5Z 9j m j k {q)q k - U{q ), (12) 

1 j,k=l 

where the scalar inertial terms m,jk ■ (S 2 ) n —>• R 1 satisfy 
the symmetry condition rrijk (q) = mkj (q) and the quadratic 
form in the time derivatives of the configuration is positive- 
definite on (S 2 ) ra . 

Corollary 1: The Euler-Lagrange equations for the La¬ 
grangian given by dT 2 l> with the generalized force fi are 


m u(q)di + (J 3 - q i q[)'^2rn ij {q)q :j + m u (q) ||g;|| 2 q t 


3 =i 


+ (h - ml) {Ad A) + 


-fi\= 0, (13) 


dU(q ) 

dqi 

for i £ {1,... ,n}, where Fi(q,q) £ R 3 is 

n i q n 

Ad,d) =^2rh ij (d)q j - -Q- dj m jk(d)dk- 

j =1 _ 1 j,k— 1 

Proof: See Appendix IBl ■ 

Note that the third term of (fl3l > represents the centripetal 
acceleration to ensure that the unit-length constraint ||f / 7 || = 1 
is always satisfied. If the inertia terms rriij are independent 
of the configuration q, then Fi(q, q) = 0 . 


B. Euler-Lagrange equations in terms of ( q , u) 

An alternate expression for the Euler-Lagrange equations 
is now obtained in terms of the angular velocities introduced 
in ©. We express the action integral in terms of the modified 
Lagrangian function 

L{d,A = L (d,S{oj)d), (14) 


where the kinematics equations are given by ©. We use the 
notation ut = (wi,..., uj n ) £ R 3 ”, and we view the modified 
Lagrangian L(q, u>) as being defined on the tangent bundle 
T(S 2 ) n . 

Let the modified action integral be © = J f J L dt. Its 
infinitesimal variation can be written as 


5© 



9L(q,A 

duji 


■ Suit + 


9L{q,ui) 

dqt 



dt. 

(15) 


Next, we derive expressions for the infinitesimal variation 
of the angular velocity vectors. Lrom ©, 


SuJi = S(Sq l )q i + S(qi)Sqi, 

Substituting <[T|» and © and rearranging. 


Sui = {SA)di) X di+ S(qi)(S(fi)di + S("/i)di)- 


Expanding each term and using the fact that qi ■ qi = 1 and 
qi • di = di ■ = 0 , this reduces to 

6u)i = (qi ■ Adi + 7* - (qi • ji)di- 

Substitute © to obtain 

SuJi = (ji ■ (S(Ui)di))di + (h - didJ)A 
= didl(S(AA + (h - didJ)A 

The matrix q,f/[ corresponds to the orthogonal projection 
along qi. But, as both of and uii are orthogonal to r/,, 
S(A ljJ i is already parallel to q,. Therefore, 

5ui =-S(Ali +(h-didj)ii- ( 16 ) 

The infinitesimal variation of w, is composed of two parts: 
the first term —S(uji)yi = 7 ^ x w, is parallel to qi, and it 
represents the variations due to the change of q.j ; the second 
term corresponds to the orthogonal projection of 7 ,; onto the 
orthogonal complement to qi, and it is due to the time rate 
change of the variation of qi. 

We now substitute © and ( fl~ 6 t into (fl5l ). and simplify the 
result to obtain the Euler-Lagrange equations expressed in 
terms of (q,ui)- 

Proposition 2: The Euler-Lagrange equations on (S 2 ) ra 
for the Lagrangian given by (fTTb with the generalized force 

fi are 


(I -did!) 


d (dL(q,A 
dt \ dcui 


2 SA) 


dL(q,A 

dtoi 


~ s (d*) + f^j = i = l,...,n. (17) 

Proof: See Appendix |C] ■ 

This form of the Euler-Lagrange equations on (S 2 ) n , 
expressed in terms of angular velocities, can be obtained 
directly from the Euler-Lagrange equations given in (fTTb by 
viewing the kinematics © as defining a change of variables 
from q to oj. This establishes the equivalence of the Euler- 
Lagrange equations in terms of (q,u>) (fTTb and the Euler- 
Lagrange equations in terms of (q, q) (fTTb . 










Next, we consider the important case that the kinetic 
energy is a quadratic form as in Substituting © 

into (ITIl i. and using the fact that ujJS( qi) T Si^qfjuji = 
cof (I 3 — qiqf) LOi = t*jJLOi, the modified Lagrangian can be 
expressed as 

1 n 

L(q,u) = -^ 2 u}fmii(q)uJi 
^ 2=1 
^ n n 

w ? S (qi) T mij(q)S(qj)ujj ~U(q). (18) 

i=i j =1 

oAi 

Substituting this into (IT7l) yields the corresponding Euler- 
Lagrange equations as follows. 

Corollary 2: The Euler-Lagrange equations for the mod¬ 
ified Lagrangian given by ( IT8l ) with the generalized force /,; 
are 

n 

mu{q)LOi + y ^S(q i ) T m ij (q)S(qj)u) : j 
f=i 

-m ij (q)S(q i ) ||wjf qj 
+ S(q i )^F i (q,uj) + ^^--fi | = 0, (19) 

for i = 1,..., n, where Fi(q , uj) G R 3 is 

n 

Fi{q,u) =^ 2 m i j{q)S(uj j )q J 
j=i 

Tor \Tor \ \ dm jk{q) 

~o 2 ^ 2 ^( q J S ^ S(uj k )q k ) . 

j=l fc=1 * 

Proof: See Appendix |C] ■ 

Similar to Corollary 1, if the inertial terms are independent 
of the configuration, then Fi(q,uj) = 0. This version of the 
Euler-Lagrange differential equations describe the dynamical 
flow (g,w) G T(S 2 ) n on the tangent bundle of (S 2 ) n . 

IV. Hamiltonian Mechanics on Two-Spheres 

We will now introduce the Legendre transformation and 
then we derive Hamilton’s equations for dynamics that 
evolve on (S 2 )'"'. The derivation is based on the phase space 
variational principle, a natural modification of Hamilton’s 
principle for Lagrangian dynamics. Two forms of Hamilton’s 
equations are obtained. One form is expressed in terms of 
momentum vectors (pi,..., p n ) € T*(S 2 )" that are con¬ 
jugate to the velocities (<ji,... ,q n ) € T g (S 2 ) n , where q G 
(S 2 )" . The other form of Hamilton’s equations are expressed 
in terms of momentum vectors (tti ,..., n n ) G K 3 ™ that are 
conjugate to the angular velocities (w 1,..., ui n ) G R 3ra . 

A. Hamilton’s equations in terms of (q, ff) 

As in the prior section, we begin with a Lagrangian 
function L : T(S 2 ) n —» R 1 , which is a real-valued function 
defined on the tangent bundle of the configuration manifold 
(S 2 )". The Legendre transformation of the Lagrangian func¬ 
tion L(q , q) leads to the Hamiltonian form of the equations 
of motion in terms of a conjugate momentum vector. Lor 


qi G S 2 , the corresponding conjugate momentum fi, lies in 
the dual space Pf S 2 . We identify the tangent space T gi S 2 
and its dual space T*. S 2 by using the usual dot product in 
R 3 . More explicitly, the Legendre transformation is given by 


Pi- qt = 


dL(q,q) 

din 


■ it, 


for any rf G R 3 orthogonal to qi. Since the component of //., 
parallel to qi has no effect on the inner product above, the 
vector representing \Xi is selected to be orthogonal to q -,,; that 
is Hi is equal to the projection of c)L q onto the tangent 
space T g; S 2 . Thus 


Pi = (h - qiqj) 


dL(q,q) 

dqi 


( 20 ) 


We assume that the Lagrangian function has the property 
that the Legendre transformation is invertible in the sense 
that the above n algebraic equations, viewed as a mapping 
from T g (S 2 )” to T*(S 2 )”, is invertible. Since these tangent 
and cotangent spaces are embedded in R 3 ", we can view 
the Legendre transformation as being the restriction of a 
mapping from R 3 ” to R 3 ” that is invertible. 

The Hamiltonian function H : T*(S 2 ) n -A R 1 is given by 


n 

H{q,p) = '^2m-qi~ L(q,q), ( 21 ) 

i=1 


where the right hand side is expressed in terms of (q, /j) 
using the Legendre transformation ( l20l >. 

The Legendre transformation can be viewed as defining 
a transformation (q, q) i-a (q, ft), which implies that the 
Euler-Lagrange equations can be written in terms of the 
transformed variables; this is effectively Hamilton’s equa¬ 
tions. However, Hamilton’s equations can also be obtained 
using Hamilton’s phase space variational principle, and this 
approach is now introduced. 

Consider the action integral in the form. 



Integrating by parts and using the fact that the variation Sq 
vanishes at to and tf, the infinitesimal variation of the action 
integral is given by 


r 

to 


»=1 

+ ( qi — 


Pi - 

dH(n,p) 

dpi 


dH(q,n) 

dqi 


■ Sqt 


■ 8 pi > dt = 0. 


( 22 ) 


Next, we derive the expression for the variations of 
According to the definition of the conjugate momenta /i, 
given by ( l20l >. we have qi ■ Hi = which implies that 
8 qi ■ Hi + Qi • dpi = 0. To impose this constraint on the 
variations explicitly, we decompose dpi into the sum of 
two orthogonal components: one component parallel to qi, 
namely Spf = qtqjSpi, and the other component orthogonal 
to qi, namely Spf 4 = (/ 3 X 3 — qiqj) 8 pi. Satisfaction of 









n 


n 


the constraint implies that qfSpf = qfSpi = —pfSqi, 
so that Spf 1 = (I 3 — qiq[) 8 ni is otherwise unconstrained. 
Substituting this and (0 into ([22l i. we obtain Hamilton’s 
equations in terms of (g, p) as follows. 

Proposition 3: Hamilton’s equations on (S 2 ) ra for the 
Hamiltonian given by (l2ll with the generalized force /', are 


Qi — (hx3 — 
Pi = —(hx3 


T .dH(q,p) 

QiQi )—a-. 

O [li 

T^ ( 9H{q,p) 


dH{q,p) 

dpt 


X {pi x qi ); 



(23) 


(24) 


for i = 1 ,... ,n. 

Proof: See Appendix IdI ■ 

Thus, equations (l23l > and (l24l > describe the Hamiltonian flow 
in terms of {q,p) £ T*(S 2 ) n on the cotangent bundle of 
(S 2 )". 

When fi = 0, any time-independent Hamiltonian is 
preserved along the solution of Hamilton’s equations, since 


dH 

dt 


dH 

~dt 


E 

i=l 


dH 

dqi 


■ qi 


dH 

dpi 


Pi 


dH 

~dt 



dH 

dpi 


x {pi x qi) 


dH 

~dt' 


Next, we consider the case where the kinetic energy is a 
quadratic function of the time derivatives of the configuration 
so that the Lagrangian is given by (ILB . The conjugate 
momentum vector is defined by the Legendre transformation 

n 

Pi = m u (q)qi + (J 3 - qiqj) Y m l] {q)q j . 

J=i 

We assume that these algebraic equations, viewed as a 
linear mapping from (gi,...,g„) £ T g (S 2 ) n C R 3n to 
(pi,..., p n ) £ T*(S 2 ) n C R 3n , can be inverted and 
expressed in the form 

n 

Qi = {h - qiqj) Y m lM O/fn ( 25 ) 

j=i 

where rn^ : (S 2 ) ra —> R 3x3 . There is no loss of gen¬ 
erality in including the indicated projection in the above 
expression since the inverse necessarily guarantees that if 
(p 1 : ...,p n ) £ T*(S 2 ) n C M 3n then (q 1: ...,q n ) £ 
T q {s 2 r c R 3n . The Hamiltonian can be expressed as 

1 n 

H{q,p) = - Y E m \ k (o)Pk + U{q). ( 26 ) 

j,k =1 


Corollary 3: Hamilton’s equations for the Hamiltonian 
given by (l26l> with the generalized force /, are (125 b and 

n 

Pi =Y ( m ij(Q^j) X X Qi) 

3 =1 



- (h - ml) - /.) . (27) 

Proof: See Appendix IdI ■ 

Hamilton’s equations (l25l > and (ITTI) describe the Hamil¬ 
tonian flow in terms of (q,p) £ T*(S 2 )" on the cotangent 
bundle of (S 2 ) n . 


B. Hamilton’s equations in terms of {q,n) 

We now present an alternate version of Hamilton’s equa¬ 
tions using the Legendre transformation of the modified La¬ 
grangian function L{q, u>) to define the conjugate momentum 
vector. The Legendre transformation (wi,..., w n ) £ R 3 " —$■ 
{n ±,..., 7 r n ) £ R 3n is defined by 


T dL(q,u) 
^ = ^- q i q i)—g—- 


(28) 


Here n^ £ R 3 is viewed as conjugate to Wj £ R 3 , i = 
1,..., n. We use the notation 7 r = (7Ti,..., n n ) £ R 3n . 
We assume that the modified Lagrangian function has the 
property that the Legendre transformation is invertible in 
the sense that the above algebraic equations, viewed as a 
mapping from R 3 ” to R 3n , is invertible. 

The modified Hamiltonian function is given by 


H{q,n) = Y' K o -ujj -L(q,ui), (29) 

3 =1 


where the right hand side is expressed in terms of (g, 7 r) 
using the Legendre transformation (f28l> . 

Consider the modified action integral of the form. 



dt. 


Take the infinitesimal variation of © and integrate by parts 
to obtain 



where we use the fact that (/ 3 X 3 — Qi/lT)^ = 7r* since 7 r* is 
orthogonal to qi by the definition (l28l >. 

The orthogonality condition ni ■ qi = 0 also implies that 
Sqi ■ ni + qi ■ Sni = 0. To impose this constraint on the vari¬ 
ations explicitly, we decompose 5ni into a component that 
is parallel to qi, namely Snf = q,(f[dn t , and a component 
that is orthogonal to qi, namely Snf 1 = (/ 3 X 3 — qiqj) 8 ni. 
From the above constraint, we have qjSni = — irJ'Sqi = 
-n?S(ii)qi = njS{q i ) r ) i . Therefore 8 nf = q t qjSn t = 
qin? S(qi)ji. 

Proposition 4: Hamilton’s equations for the modified 
Hamiltonian given by ( [29b with the generalized force f are 


Qi = ~S{Qi) 


9H{g,n) 

dni 


(31) 










TTi 


~S(qi) 


dH(q,ir) 

dqi 


dH(q,n) 

dni 


x 7 u + S(qi)fi, 


(32) 


for i = 1 ,... ,n. 

Proof: See Appendix [E] ■ 

Thus equations (f3Tb and (f32l) describe the Hamiltonian flow 
in terms of (q,7r) £ T*(S 2 ) n on the cotangent bundle of 
(S 2 )". 

When fi = 0, any time-independent modified Hamiltonian 
is preserved along the solution of Hamilton’s equations, since 


dH 

dt 


dH 22 

dH 

dt + dqi 

1=1 

’ Qi + t;— ' 

U7Ti 

dH 2A 

( dH 

~dt + ^d7 z 

2=1 

' X7F 


dH 

~dt‘ 


We now consider the important case where the kinetic 
energy is a quadratic function of the angular velocities in 
the form that arises from the Lagrangian given by (IT4l >. The 
conjugate momentum is defined by the Legendre transfor¬ 
mation 


TTi =mn(q)wi + 22 S(q i ) T m ij (q)S(q j )u>j. (33) 

i=i 

i#* 


We assume these algebraic equations, viewed as a linear 
mapping from (wi,... ,ui n ) £ M 3 ™ to (7Ti,... ,n n ) £ M 3ra 
can be inverted and expressed in the form 

n 

uj i = '^2 m ij(q) n P ( 34 ) 

j=i 

where mfj : (S 2 ) n —» IR 3x3 . The modified Hamiltonian 
function can be expressed as 


ii it 

H{q, 7t) = - 22 22 + U A)- (35) 

*=1 3 = 1 


Corollary 4: Hamilton’s equations for the modified 


Hamiltonian given by (l35l> are 



Qi = ~S(qi) | 

^2 m ij(9) 7L j, 


(36) 

i = -S(q») | 

\ i a T , , . 

1 L ”i + 

{ j,k =1 

dU(q) 

dqi 

! 

+ < 22 m ij( q ) 7r j f x Tti + S(qi)fi, 


(37) 

for i = 1,... ,n. 

Proof: See Appendix [E] 


■ 


V. Dynamics on Chain Pendulum 


A chain pendulum is a connection of n rigid links, that are 
serially connected by two degree-of-freedom spherical joints. 
We assume that each link of the chain pendulum is a rigid rod 
with mass concentrated at the outboard end of the link. One 
end of the chain pendulum is connected to a spherical joint 


that is supported by a fixed base. A constant gravitational 
acceleration acts on each link of the chain pendulum. This 
may represent a spherical pendulum (n = 1), or a double 
spherical pendulum (n = 2) as special cases. 

We demonstrate that globally valid Euler-Lagrange equa¬ 
tions can be developed for the chain pendulum, and they 
can be expressed in a compact form. The results provide an 
intrinsic and unified framework to study the dynamics of a 
chain pendulum system, that is applicable for an arbitrary 
number of links, and globally valid for any configuration of 
the links. 

The mass of the <-th link is denoted by rn t and the link 
length is denoted by U. For simplicity, we assume that the 
mass of each link is concentrated at the outboard end of 
the link. An inertial frame is chosen such that the first two 
axes are horizontal and the third axis is vertical. The origin 
of the inertial frame is located at the fixed spherical joint. 
Each of the chain links has a body-fixed frame with the 
third axis of the body-fixed frame aligned with the axial 
direction of the link. The vector q\ £ S 2 represents the 
direction from the fixed base to the mass element of the 
first link, and the vector q t £ S 2 represents the direction 
from the (i — l)-th spherical joint to the concentrated mass 
element of the v'-th link. Thus, the configuration of the chain 
pendulum is the ordered n-tuple of configurations of each 
link q = (qi,...,q n ) £ (S 2 ) n , so that the configuration 
manifold is (S 2 ) n . The chain pendulum has 2 n degrees of 
freedom. 

Let Xi £ R 3 be the position of the outboard end of 
the i-th link in the inertial frame; it can be written as 
x i = ljq : j. The total kinetic energy is composed of 

the kinetic energy of each mass: 

1 n i 

T(q,q) = y:J2 m ^J2 l Aj\\ 2 - 

i=1 3= 1 

This can be rewritten as 
1 n 

T(q,ui) = — 22 MijliljUjJS(qi) T M i jl i ljS(q j )u}j, (38) 
i,j =i 

where the real inertia constants are given by 

Mij = I 22 mk I ’ *>3 = ''' ’ n - 

\fe=max{i,j} / 

The potential energy consists of the gravitational potential 
energy of all mass elements. The potential energy can be 
written as 

n n n 

u{q) =22 mi 9 e 3 x i =22'^2 m A l i e 3qi- (39) 

i= 1 i—1 j=i 

The modified Lagrangian function L : T(S 2 )" — > R 1 of the 
chain pendulum is given by L(q, to) = T(q,co) — U(q) from 
([31 and I®. 

Also, suppose that there exist a control torque r £ K 3 
acting on the spherical joint connecting the fixed base and 
the first link, and a disturbance force d £ M 3 acting on the 





tip of the last link. The corresponding virtual work is given 
by 

n n 

6 W = 71 • t + Y hS{li)Qi ■ d = 71 ■ r + • kS{qi)d. 

2=1 2=1 

Therefore, the generalized forces are given by /i = r + 
hS(qi)d, and f 3 = l 3 S{q 3 )d for j > 2. 

Substituting this into CD . the Euler-Lagrange equations 
for a chain pendulum are given by 

n 

MuliUi + Y', (qi)S(qj)ibj 

j =i 

n n 

- ym il 2 S{ qi ) qj ~2mjgliS(e 3 )qi = S(qi)f x . 


j =i 


(40) 


Similarly, the Legendre transformation is given by ( |3P . 
and from (l36b . (f3Tb . Hamilton’s equations can be written as 


qi = -S(qi) ^ ) > 

1 d 


(41) 


TTi = -S{qi) < Y n J M jk(qW + Y rn o9h e 3 


2 

+ X! ^ X Vi + S (qi)fi- 

,i=! 


(42) 


These are remarkably compact considering the complexity 
of the dynamics, and they are well structured compared with 
the equations of motion expressed in terms of angles. 

This mathematical model may be applied to a wide class of 
other dynamical systems, such as articulated robotic systems. 
As they are developed for an arbitrary number of links, 
they are readily extended to finite-element approximations 
of cables or slender rods after augmenting the potential with 
an elastic potential term. The proposed global formulations 
avoid singularities associated with local coordinates. 

Appendix 

A. Hat map 

Several properties of the hat map are summarized as 
follows. 

S(x)y = x x y = —y x x = —S(y)x, (43) 

S(x) 2 = -{x T x)I 3x3 + xx T , (44) 

S(x) 3 = —(x T x)S(x), (45) 

x ■ S(y)z = y ■ S(z)x = z ■ S(x)y, (46) 

S(x)S(y)z = (x ■ z)y - (x ■ y)z = ( yx T - x T yI 3x3 )z, 

(47) 

S(x xy) = S(x)S(y) - S(y)S(x) = yx T - xy T , (48) 

for any x,y,z G R 3 . 


B. Proof of Proposition 1 

Substituting £7) and © into £7}, and rearranging with ( I46| i. 

9L(q,q) 


S& = J E ’ ( S (9i) 


dtp 


+ .*.(S( 4 )®^ + S(«)“M'| Ut 


dqi 


dqi 


Integrating the first term on the right by parts, and using the 
fact that the variation vanishes at to and f/, it reduces to 

XlK v- f tf ( 0 , ^\ d dL{q,q ) dL{q,q)\\ 

se= YX n sfe) — &r)) dt 

The virtual work is written as 

n n 

6 W = Yh- S( di)qi = Y" 1i ' s ^i)fi- 

■i =1 i =1 

According to Lagrange-d’Alembert principle, 6(3 = —6W 
for all continuous variations 7 , : [to,t/| —>• R 3 , that satisfy 
( 7 * ■ qi) = 0. The fundamental lemma of the calculus of 
variations implies that the expression — jpf — fi is 

parallel to %, or equivalently, CD- 

Next, consider the case where the kinetic energy is given 
as a quadratic form as (ITD . Substituting (THU) into CD- 

{ n n 

Y + X! "Hj X)<ij 

j=i 3 =1 


1 d , T dU{q) 


2 % EE3j m jk {q)qk 
j= 1 fc= 1 


dqi 

„T • \ _ („T-a 


-/, =0- (49) 


Since qf % = 0, it follows that ^ (qfqf) = (qf qf) + ||% || = 


T • 

qfqi 

0 ; thus we obtain 

(h - qiqf)q% = q% - (%% T )% = Qi + 11%II 2 q%- 

Substituting this into d49l ) yields CD- 

C. Proof of Proposition 2 

Substituting 0 and CD into CD and using the integration 
by parts, 


V" f tf J ^ (it T\dL(q 1 ix) 


dun 


+ s( „)“M + «^ d U t. 


dqi 


duii 


According to Lagrange-d’Alembert principle, 63 = —<5W 
for all differentiable functions 7 * : [to,tf] —> R 3 that satisfy 
( 7 i • qf) = 0 and vanish at to, and tf. This implies that the 
expression in the braces of the equation below is parallel to 
qi, or equivalently 




dL(q,u) 

duii 













+ S( qi ) + f^j } - °’ ( 5 °) 

where we used —S(qi) 2 = I 3X 3 — Qi.qJ obtained by (l44t . 

To further simplify these expressions, we derive a few 
identities. From 07} and w t ■ r/,; = 0, 

S(qi)S(uJi)S(qi) = S{qi){-uJ ql 3 + qnvj} = 0. 

From this, (0, and 08} . it follows that 

S(qi)S(qi) = S(qi){S(ui)S(qi) - S(qi)S(ui)} 

= -SiqtfS^), 

and similarly, S(qi)S(qi) = S(oJi)S(qi) 2 . Consequently, 
these results can be used to obtain 

S(qi) 2 {S(qi)S(qi) + S(qi)S(qi)} = -S(qi) 4 S(ui). 

Substituting these into 00} and using ( 1441 and ((43} repeat¬ 
edly yield ( fTTb . 

Next, the Euler-Lagrange equations for the case that the 
kinetic energy is given as a quadratic form, namely (IT9t can 
be obtained by either substituting (fl~ 8 b into (1 17 k or rewriting 
(fl3l> in terms of the angular velocity as a change of variables. 
Here, we follow the latter approach as the corresponding 
proof is more concise. 

From the kinematics equation 0 and 07}. we have 

qi = Ui X qi + u>i X (wj X qf) = -S(qi)dii - ||wi|| J <M. 


We now invoke Hamilton’s phase space variational principle 
that <5© = — 5W for all possible functions 7 , : [to,tf] —> K 3 
satisfying 7 i • qi = 0 and Spf 1 : [to,tf] —> R 3 that are 
always orthogonal to qi for i = 1,... ,n. According to the 
fundamental lemma of the calculus of variations, 


S(qi) 



dH(q,p) 

dqi 


T 

- HiQi 


dH(q,p) 

dpi 



= 0 . 


We multiply this by S(qi) and use a matrix identity to obtain 


(^3x3 



dH(q,p) 

dqi 


T 

PiQi 


dH(q,p) 

dpi 



= 0 . 


Since both terms multiplying Spf 1 in the above variational 
expression are necessarily orthogonal to q % , it yields (03} . 

We now determine an expression for /i,;. The above 
equation only determines the component of /i, that is normal 
to qi. The other component of /i, that is parallel to q, is 
derived by taking the time derivative of qi ■ /j,i = 0 to obtain 
qi ■ fii = —qi ■ Hi. Thus, Hi is obtained by the sum of two 
components as 

n T\ ( dH T dH\ T 

Pi = ihx3 - QiQi ) 1—^7 + PiQi J - (Mi ViM, 

which is reduced to 04} via (03} and 07}. 

Next. 07} can be derived by substituting (06} into (04} 
directly and rearrange it. 


Substituting this into 03}, and rearranging it with I 3 x 3 

ml = -S(qi) 2 , 

n 

-S(qi)^mu(q)ui + S(q i )'^2m ij (q)(-S(q j )dj j - \\uj\\ 2 qj ) 

7=1 

+ S(qi)(Fi{q,u) + dU ^ - /i)} = 0 . 


E. Proof of Proposition 4 

Using Sni = dn^+Sn^, Snf = q t qf Sni = q^JS(q i )') i , 

and qjUi = 0 , the variation of the action integral given at 
00 } can be rewritten as 

<5© = Y ~][ { - 7fj + S(ui)iTi - S(qi ) ^ (51) 


In the above equation, the left hand side becomes zero when 
the expressions in the braces are either zero or parallel to 
qi. However, the second case is not possible as they are 
perpendicular to qi by the definition. This yields 09} . 


D. Proof of Proposition 3 

Substituting Spi = Spf + Sp 44 and 0 into 02}, and 
rearranging it with qfSpf = — pfdqi , 




pf 


't 0 


+ q^ - 


dH(q,p) 

dH(q, p) 

dpi 


■ Sqt 


+ {hx3 - ml) ii - 


Sp? 

dH(q,p) 

dpt 


= £/ 

i=i 


S(qi) -Mi - 


dH(q,p) 

dqi 


■ Pi Q: 


■ tPi j dt, 

t d H(q,p) 
dpi 


. I ■ ,r T\dH(q, p) 

+ ( Qi — (13x 3 — QiQi )—- 


6 p? } dt. 


S(qi)m 


T dH(q, 7r) 


diTi 

Ui - (J 3 x 3 - ml) 




T \dH(q,n) 


diri 


■Snfdt. (52) 


According to Hamilton’s phase space variational principle, 
<5© = —6W, the expression at the first pair of braces of the 
above equation, added with S(qi)fi, should be parallel to q t , 
or equivalently, 

(hx3 - ml){ - Tti + S(uJi)TTi - 

+S(qi)7nql dH ^ ?r) + S( qi )fi} = 0. (53) 

Also, the expression at the second pair of braces of 02} is 
already parallel to qi to yield 


Ui — {I 3 x3 


T\ 

mi ) 


dH{q, ?r) 
dni 


-S(q t ) 2 


dHjq, tt ) 

dni 


Substituting the second equation into 0 yields 01} . 





















Using the facts that (J 3x3 - qiqf)S{qi) = -S(qi ) 3 = 
S{qi) and (I 3x3 -qiqJ)S(ui)-Ki = -£(&){& X (w* X7T*)} = 
-S(qi){(qi ■ 7 Ti)u>i - (qi • = 0. ((53]) reduces to 



However, this is incomplete since it only determines the 
component of j ij that is perpendicular to qi. The component 
of TTi that is parallel to qi is determined by taking the time 
derivative of qi -TTi =0 to obtain qt --hi = —qi-iTi. Therefore, 
qi,qj 71 i = —qinfqi. By combining these. 



But, —S(qi) 2 TTi = ni since 7r* is normal to qi. This yields 



Next, (l36l > is obtained by substituting the angular velocity 
(ful l into the kinematics equation <[3}, and ( 13 71 is derived by 
substituting the Hamiltonian (l35l > into (l32l >. 
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